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Abstract 

Monte Carlo simulations were performed to investigate the possibility to add a 
DIRC detector to the WASA detector at the COSY accelerator. A statistical 
method for pattern recognition is presented and the possible angle resolution and 
velocity precision achieved are discussed. 

1 Introduction 

It is planned to measure rare decays of rf mesons with the WASA at COSY 
detector [1], [2]. The foreseen production reaction is 

p + p^r/' + p + p (1) 

where the two protons will be measured and the rj' will be identified via 
missing mass technique. The cross section for this process is small, because 
no nucleon resonance couples strongly to the rj'p channel, in contrary to r] 
production. The rf peak will be on top of a huge multi-pion background. 
The only chance to identify the 77' is to measure the four momentum vectors 
of the two protons with high accuracy. Although their emission angles 
(9, 4>) are measured with sufficient precision, their kinetic energies are 
poorly determined by the present setup. A possible option to achieve 
sufficient resolutions seems to be a DIRC detector (Detection [of] Internally 
Reflected Cherenkov [light]) [3] in the forward detection system of WASA. 
Cherenkov detectors are widely used in high energy physics. They can have 
different applications: particle identification, energy measuring, threshold 
detector for detecting specific type of particles in a definite energy range. 
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Our aim is to measure the velocity of protons in the kinetic energy range 
from T = 400 MeV to T = 1000 MeV with a precision better than 

= 0.5% 

2 DIRC Design 

With the Monte Carlo package GEANT3 [4] different radiator materials 
were studied: LiF, CaF 2 , quartz (fused silica) and plexiglas. All materials 
have different refractive indices which are smooth functions of the wave 
length in the visible range (see Fig. [1]). Quartz has the largest values for 
the refraction index n at any given wavelength A compared to the two other 
materials, n varies from 1.60 in the UV range to 1.45 in the IR range. The 
transmittances are rather independent from the wavelength above a certain 
"threshold" wavelength as is shown in Fig. [2] (data taken from Ref. [5] 
and [6]). The radiator should have a large refractive index and high 
transmittance to measure particles with "low" kinetic energy. LiF has 
higher transmittance in the UV range compared to CaF 2 and quartz. This 
gives the possibility to build a Cherenkov detector with LiF as radiator and 
a proper photocathode for detecting photons in the UV range. However, in 

the visible range, where the transmittance of LiF and quartz are 
compatible, quartz has the larger refractive index. We therefore omit CaF 2 

and LiF in the further studies. 
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Figure 1: The refractive index n as function of the wavelength for the indi- 
cated materials (from Refs. [5,6]). 
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Figure 2: Radiator transmittances (solid curves) [5] and quantum efficiencies 
of photocathodes (dotted curves) [7-10] as function of the wave length. The 
bulk transmittances are given for 10 mm thickness and reflection losses at 
the surface result in an additional reduction of the transmittance by 10%. 



The quantum efficiencies for different photocathodes are also shown in Fig. 
|2] (from Refs. [7-10]). Csl and CsTe are used for photon detecting in the 
UV range. Bi-alkali (KCsSb) and multi-alkali (GaAsP) photocathodes are 

used for the near UV and the visible light range. These types of 
photocathodes have larger and broader distributions of the corresponding 
quantum efficiency function. This could compensate a lower number of 
emitted photons in this range of wavelengths. In a first step, simulations 
were performed for a RICH type Cherenkov detector with different 
radiators and photon detectors. However, a DIRC type detector was found 
to be favorable due to following reasons: energy losses in the detector, the 
broad range of velocities to be measured in the proximity to the Cherenkov 
threshold and the availability of very fast photon detectors. Based on the 

above discussed properties of the materials we limit ourselves in the 
simulations of DIRC detectors to fused silica and UV transparent plexiglas 

as radiator and photomultipliers with bialkali photocathodes. 
We will now concentrate on the DIRC detector. The design differs from the 

Babar DIRC [3,11]. It consists of radiator bars with dimensions 2x6x50 
cm 3 and a focusing element with a cylindrical mirror (Figs. [3] and H]) . The 
dimensions were for fused silica whereas for plexiglas an increased thickness 
of 4 cm was assumed. The individual modules are arranged in two planes, 
one above and one below the beam axis. The normal axis to each plane is 
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tilted by 20° away from the beam axis (see Fig. [3]). Emitted Cherenkov 
photons create an image on a pixelised screen. Each pixel is a bi-alkali 
photocathode with dimension 5x5 mm 2 . 




Figure 3: Schematic side view of the detector. The radiator bars are inclined 
by 70 degree with respect to the beam axis. A simulated proton is indicated 
with an arrow, the Cherenkov light as the zigzag and arrow. The focussing 
bar (mirror) is on top. 



3 Monte Carlo Results 

A single event image consists of two arcs, which are the results of the 
internal reflection from the left and right side of the radiator bar. The 
shape of the arcs depend on the kinetic energy and particle species. The 
relative positions of the arcs depend on the horizontal angle of the detected 
particle track: the larger the horizontal angle, the larger the distance 

between arc centers [3]. 
We have simulated Cherenkov images for protons with 600 MeV and 1000 
MeV kinetic energy. We applied a left handed coordinate system with the 
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Figure 4: Assembly of several bars to one module. Two such modules are 
foreseen, one above and one below the beam. 



z-axis in the beam direction, x — z defines the horizontal plane and y is 
vertical to this plane. The angle between the ejectile momentum and the 
(x-z) plane is # vert and the angle between the momentum projected in the 
plane and the beam direction is #horiz- Thus 6*horiz = arctan(za;/p z ) and 
#vert — arcsin(p y /p). Simulations were performed for track directions 
#horiz ~ 8.5° fixed and varied and vice versa: 9 VCTt ~ 8.5° fixed and 6*horiz 
varied. They are shown in Fig. [5] 




Figure 5: Left: Scatter plot of the position of detected photons from 100 
simulated events for protons with kinetic energy T— 600 MeV. The abscissa 
and ordinate are the horizontal and vertical position numbers on the photo- 
cathode. Right: Same as left but for T=1000 MeV 

The average number of emitted photons depends on the path through the 
radiator end the kinetic energy. Fig. [6] shows the dependence of the number 
of detected photons on the horizontal track direction for protons with 
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kinetic energy T = 600 MeV. The number of detected photons drops 
sharply when the particle passes through the edges of the radiator bars. 
This effect results from the detector construction. 
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Figure 6: The average number of the detected photons as function of the 
horizontal angle for protons with T = 600 MeV. 



Furthermore, image distortion is also observed. Smeared arcs or even four 
arcs could be observed when the particle goes through the edges of two 

neighboring bars. 

4 Event Reconstruction 

The following method is used for event reconstruction. A collection of 
events with different particle species as well as energies and directions is 
simulated. Then a second set of simulated data named detected are 
"compared" to the simulated data. The comparison is realized by the 

following function F ([5} 

F(k,(E,P)) = J2det(i,j)-sim(i,j). (2) 

Here k denotes the particle species, E, P the energy and momentum of the 
particle, det(i,j) is the number of detected photons in (i, j)-th 
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photoelement in the experiment and sim(i,j) the simulated average 
number of detected photons in (i,j)-th photoelement. The expression is a 

maximum when the experimental image on the screen is similar to the 
simulated one. This method provides an opportunity to determine (3 and 
the direction of a particle if the particle type is known. The particle species 
information has to come from other sources. 
The possibility of track reconstruction for protons with T=600 MeV is 
demonstrated in Figs. [7] and [HJ The particle direction is extracted under 
the assumption that the particle species and its energy are known. 
Horizontal and vertical angles relative to the beam direction are 
reconstructed given that one of both is fixed. The small widths (a) of 
Gaussians fitted to the data of ~0.3° and ~0. 5° indicate the sensitivity of 
the method. Similar results were obtained for protons with T=1000 MeV. 



500 



400- 



Q_ 

llT 



300- 



200- 



100 





- w 8 - 531 '; 

: 




i 


; J 


I • 




< k — 


2 4 6 8 


10 12 14 16 



500- 



400- 



Q_ 

llT 



300 



200- 



100 



=8.531 



V 



.( ) 



4 8 12 16 



Figure 7: Track reconstruction. Left: F as function of the vertical angle and 
for fixed horizontal angle. The results (dots) are connected by lines to guide 
the eye. The smooth curve is a fit with a Gaussian plus a constant. Right: 
Same as left but horizontal and vertical angles are interchanged. 



Fig. [8] shows the /^-resolution for protons. Here we have assumed that the 
angle measurement is performed with other detectors. The method gives a 
resolution (crp/(3) for protons better than 0.3%. The results obtained can 
be transformed into the missing mass resolution. For an incident proton 
beam with 3.3 GeV/c momentum a missing mass resolution of a mm ~ 0.5 
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Figure 8: Relative /3-resolution for protons as function of (3. The results are 
shown by dots, the solid line is a linear fit to the dots. 



MeV is achieved. 

5 Summary 

We have designed a DIRC detector to measure protons with kinetic energy 

from 400 to 1000 MeV. Monte Carlo simulations for such a device have 
been performed. A method is discussed to deduce four momentum vectors 
of the protons. The method gives reliable information for the event 
reconstruction: velocity-resolution for protons is better than 0.3%. The 
routine is able to reconstruct tracks and identify particles when additional 
information is provided by complementary detectors. The large data base 
of simulated events is a disadvantage of this method; however, optimization 
of the algorithm will increase the speed. 
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